Gap and screening in Raman scattering of a Bose condensed gas 
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We propose different spectroscopic methods to explore the nature of the thermal excitations of a 
trapped Bose condensed gas: 1) a four photon process to probe the uniform region in the trap center: 
2) a stimulated Raman process in order to analyze the influence of a momentum transfer in the 
resulting scattered atom momentum distribution. We apply these methods to address specifically 
the energy spectrum and the scattering amplitude of these excitations in a transition between 
two hyperfine levels of the gas atoms. In particular, we exemplify the potential offered by these 
proposed techniques by contrasting the spectrum expected, from the non conserving Bogoliubov 
approximation valid for weak depletion, to the spectrum of the finite temperature extensions like 
the conserving generalized random phase approximation (GRPA). Both predict the existence of the 
Bogoliubov collective excitations but the GRPA approximation distinguishes them from the single 
atom excitations with a gapped and parabolic dispersion relation and accounts for the dynamical 
screening of any external perturbation applied to the gas. We propose two feasible experiments, one 
concerns the observation of the gap associated to this second branch of excitations and the other 
deals with this screening effect. 

PACS numbers: 03.75.Hh,03.75.Kk,05.30.-d 



INTRODUCTION 

Nature of the elementary excitations 

The experimental discovery of the condensation of a 
Bose gas has confirmed the existence of the phonon-like 
nature of the collective excitations 0, Q . The obtained 
measured energy spectrum not only is gapless as stated 
from the Hugenholtz-Pines theorem but is also in per- 
fect agreement with the prediction of the Bogoliubov ap- 
proach at zero temperature [3||. However, a second fun- 
damental question arises as to whether these collective 
excitations are the elementary building constituents for 
the normal part of the fluid as assumed in the Bogoli- 
ubov approximation. Most standard textbooks rely on 
this quasiparticle hypothesis in order to determine the fi- 
nite temperature gas properties [3, HI- In contrast, in the 
theoretical description of a plasma, distinction is made 
between the elementary excitations (ions) and the col- 
lective ones (plasmons). As discussed in previous works 
there are no fundamental reasons to exclude this 
distinction also in a Bose gas. 

Precisely, suppose a bulk gas of total and condensed 
densities n and no embedded in an volume V where 
atoms of mass m interact through the s wave channel 
with a scattering length a. The Bogoliubov approxi- 
mation predicts that the elementary excitations of mo- 
mentum k are phonon-like with a dispersion relation 
given by ef k = ^2gn e k + e 2 . where e k 



and g = 4irah 2 /m. Nevertheless, its non conserving 
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property (violation of the mass conservation law) [10( 
restricts its validity for a weakly depleted Bose gas 
and thus limits its use to low temperature. As op- 
posed to that, the so-called generalized random phase ap- 



proximation (GRPA) or equivalently the time-dependent 
Hartree-Fock (TDHF) approximation is instead conserv- 
ing and valid for the whole range of temperature. This 
alternative approach distinguishes explicitly these collec- 
tive phonon-like excitations from the atom-like elemen- 
tary excitations with the rjarabolic dispersion relation 

IT III . The constant term 



if = £k + ff(2n - n <5 k ,o) 



corresponds to the Hartree and Fock (HF) mean field en- 
ergy part and takes into account the absence of exchange 
interaction energy between condensed atoms. Therefore, 
an energy gap exists between the thermal and condensed 



atoms efjf - ef „ = 9 n o + «k 



HF 



Superfluidity due to total screening 

Another important reason to discriminate among the 
various theoretical approaches is to have an improved 
understanding of the superfluidity phenomenon. More 
precisely, we would like to answer the following question: 
Why, from a kinetic point of view, a superfluid can re- 
main in a metastable motion without converting its ki- 
netic energy into heat? Many explanations have been 
provided but, according to [J], the situation is not en- 
tirely clear as far as kinetic theory is concerned. 

Instead, the equilibrium aspects based on the ensemble 
approach of the superfluid phenomenon of a Bose con- 
densed gas are well understood. Using the « ensemble 
which breaks the U(l) symmetry associated to the par- 
ticle number conservation, one can describe the super- 
fluid motion (condensed mode) relatively to the normal 
fluid (non condensed modes) [9|. Such a relative motion 
should not be considered as an equilibrium state but as a 
metastable state possibly subject to relaxation of a state 
of lower energy. Unfortunately, the ensemble approach 
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does not explain the physical reasons for such a metasta- 
bility. It just tells that an artificial breaking of symmetry 
allows you such a description. Only a non-equilibrium 
treatment can provide these explanations and therefore 
confirm the validity of the assumptions used in the en- 
semble approach. 

The kinetic theory so far developed in the Bogoliubov 
approximation allows for such a metastability in the weak 
depletion limit [14]. Particle exchange between the nor- 
mal and superfluid are regulated through a balanced Be- 
liaev process of transforming one collective excitations 
to two collective excitations. A complete different sce- 
nario appears in the GRPA as it accounts for the dy- 
namical screening of any external time-dependent poten- 
tial that perturbs the gas atoms 6, 8]. The ability of the 
macroscopic condensed wave function to deform locally 
its profile allows for a screening of any external perturba- 
tion that affects the energy transition probability of any 
atom-like excitation. In particular, under some stability 
conditions [8j , a total screening forbids individual energy 
transitions involving a condensed atom. In this sense, 
the condensed atoms are gregarious since they respond 
only collectively to a perturbation via the creation of a 
phonon-like excitation. If the external potential origi- 
nates from the presence of another thermal atom, this 
total screening prevents the binary collision between this 
thermal atom and any condensed one. Therefore, con- 
trary to the Bogoliubov approach, the metastability of 
the relative motion between the normal and super fluids 
in GRPA is explained from the absence of this exchange 
collision process. 

Nevertheless, atom exchanges between the normal and 
the super fluids should always exist in any kinetic descrip- 
tion, in particular to guarantee the process of condensate 
formation. This is the case for the GRPA, but provided 
that instability conditions are satisfied 6]. For example, 
when the relative velocity between the two fluids exceeds 
the critical velocity given by the Landau criterion, the 
total screening phenomenon disappears and the binary 
collisions become again possible. 



Experimental difficulties 

Both gap and total screening phenomena have been 
predicted to appear in a Raman transition process be- 
tween two hyperfine levels of a Kb gas, but only in the 
bulk case [6]. In comparison to other methods like ra- 
dio frequency (RF) or Bragg spectroscopy, the possibility 
of momentum transfer and the distinction between scat- 
tered and unscattered atoms enable these observations. 
However, an experimental realization is still not simple 
in the real case of a trap since the gas inhomogeneity, 
combined with the short duration of the applied coupling 
potential, leads to additional broadenings of the spectral 
lines that prevent the resolution of the gap and screen- 



ing structure. In this context, a RF spectroscopy would 
have probed the whole gas which includes thermal atoms 
of the outer and inner condensate regions. Therefore, the 
distinction between various theoretical approaches is ex- 
tremely difficult as long as the transition amplitude and 
the dispersion relation of thermal atoms have a strong 
spatial dependence. 



Setup proposals 

In this letter, we propose different methods to probe 
the atoms more efficiently than the RF spectroscopy: 1) 
the Raman scattering is a two-photon process that offers 
also the possibility to transfer the momentum q to the 
scattered atoms and observe their resulting momentum 
distribution after expanding the gas; 2) a four photon 
scattering process, where two sets of two beams cross in 
the trap center, addresses selectively the homogeneous 
region of the gas (see Fig [3]). 

We apply these methods for the case of a finite tem- 
perature trapped Bose gas in the GRPA, in a bid to chal- 
lenge the Bogoliubov approach. To this end, we propose 
two concrete experimental setups that overcome the dif- 
ficulties associated with the trap: 1) The gap is observed 
from the four-photon process; 2) The total screening is 
determined in a Raman scattering. Previous theoretical 
works argue in favor of the conserving GRPA. Nev- 
ertheless, a comparison with the non conserving Bogoli- 
ubov approximation is of relevance as long as the second 
branch of individual excitations has not been observed. 



RAMAN SCATTERING 



The GRPA approach 



In a Raman transition, we start from atoms initially 
in the hyperfine level |1) = \F — l,mp = —1). Each 
mode k is characterized by its initial population Nq and 



k#0 



l/(exp[ 



-l,k 



(j,)] — 1) and its initial plane 



wave function ip\ t \s. = exp[i(k.r — €^£t)]/y/V with the 
inverse temperature j3 — 1/fcsT and the chemical poten- 
tial \x = g(2n — no). The application of a perturbation 
coupling potential Vq(r, t) — Vr exp[i(q.r — ut)] at t > 
transfers a small fraction of them into the second level 
|2) = \F = 2,mp = 1) of internal frequency ojq. The 
determination of the second spinor component of the as- 
sociated wavefunction - 02.k(r,i) of the mode k evolves 
according to the time-dependant Hartree-Fock equation 
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FIG. 1: Diagrammatic representation of the scattering of an 
atom by an external potential. An atom of momentum k 
is scattered into a state of momentum k + q by means of an 
external interaction mediated by a virtual collective excitation 
of momentum q. 
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where we define the intercomponent coupling g 12 
47rfi 2 a 12 /m. The solution is [6J: 



^2,k(r,t) 



, L</ / ° o dt / e i (" , + i0 )( t, - t )y q (r,^)^ 1 , k (r,t) 



27ri/Ci2(q, u)'){fvui' + iO — fiwk,q) 



where fiw k ,q = £ 2 k 



MF 
-l,k 



and 



-2,k+q 



.912 {n — noS^o) is the atom mean field energy in the sec- 
ond level without the exchange term. These formulae 
resemble the one obtained from the non interacting Bose 
gas except for the HF mean field terms and the screening 
factor: 



/Ci 2 (q,w) = 1 



512 \ - N k 

V ^ hw + iO - hujk c 

k ' 



(3) 



Eq.@ is interpreted in FigfT] in terms of propagators 
whose poles determine the resonance frequencies. One 
pole is associated to the individual transition between 
atoms: uj = w k>q and the other is the zero of the screen- 
ing factor and corresponds to the collective excitations 
associated to the gas rotation in the spin space: 6ui — 
uj — luq ~ [e q — (g — gi2)n]/h for 512 ~ g. Total screen- 
ing corresponds to the singularity /Ci2(q, wo,q) —> 00 and 
prevents any single condensed atom scattering 0. 

In a bulk gas, the transferred atom density for each 
mode is obtained from n 2 ^ + q(t) = [tjj^(r, ^)| 2 ^k so that 
we deduce the total atom density[2j, 



n 2 



E 



dio 



,4sin 2 (w't/2) 



|^lVi' 2 (q^-^')(4) 



expressed in terms of the imaginary part of the 
intercomponent susceptibility function Xi 2 (q, uj) = 
l/(5i2£i 2 (q,u;)). 



The Bogoliubov approach 

These results can be compared to the one obtained 
from the Bogoliubov non conserving approximation de- 
veloped in [y, @, [H| which is valid only for a weakly 
depleted condensate. This approach implicitly assumes 
that the elementary excitations are the collective ones 
forming a basis of quantum orthogonal states for the 
description of the normal fluid. Consequently, this for- 
malism predicts no gap and no screening. The creation- 
annihilation operators c ik (t), Cj lk (t) describing the var- 
ious components in the momentum space evolve ac- 
cording to ci )k (f) = e-^v/TVo^o + u+, k e" ie ^*6i, k + 
u- tk e ie ?^b\ _ k ) and c 2 , k (*) = e- l(t * +e ^ )t c 2 , k . In this 
expression, besides the collective excitation modes of 
phonon of energy ef k , a second collective mode of ro- 
tation appears with energy e^ k = e k + (.912 — g)no- 
fx = gno is the chemical potential, 61 k is the annihi- 
lation operator associated to the quasi-particle such that 
( &t i,kVk) = l/(expG8e£ k ) - 1) and u± jk = ±((e k + 
gno)/2ef k ± 1/2) 1 / 2 . Reexpressing the intercomponent 
susceptibility 



iV 



Xia(q,w) = ^/ dte^ + ^([pf(0),p^(t)}) (5) 



in terms of the autocorrelation function of the excitation 

operator p^(t) = J2k c a,k.( t ) c p,k+q.( t )/V' we calculate 
in the Bogoliubov approximation: 



xf 2 (q> w ) 



_ <5 k ,o^Vo/2 ± U 2 ±ik /(exp(±/3ef k ) - 1) 



±,k 



V(H8u - 



"2,k±q 



(6) 



In contrast to the GRPA, Eq.© describes a spin rota- 
tion transition of the condensed fraction, one transition 
involves the excitation transfer from a phonon mode into 
a rotation mode and another the excitation creation in 
the two modes simultaneously. 



Extension to the trap 

These formulae can be easily extended to the case of 
a harmonic trap Vff(r) = X)i ma; l r i/2 of frequency u>i 
by considering the local density approximation (LDA) 
5]. For a weakly inhomogeneous gas, the population in 
each mode becomes a local quantity iV k — > A r k (r). By 
making this replacement, the thermal density ny(r) = 
Ek^o^W/V, the energies e^f (r), e^ k (r), the screen- 
ing factor /Cia(r, q, uj), the potential amplitude Vr(t) and 
71 2,k(r, t) become local quantities as well. The zero mode 
density no(r) = |^o(r)| 2 is determined from: 



fi 2 Vfoo(r) 
2m$o(r) 



T/ ff (r)+ 5 (|*o(r)| 2 + 2n T (r))= M (7) 
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FIG. 3: Selective four lasers interaction with atoms within the 
trap center region. Two lasers a + polarized along the z axis 
(1) drive the atoms through the intermediate states mj- = 
and mj? = 1 while two others along the y axis (2) drive the 
atoms within the sublevel tuf — 1. 



FIG. 2: Density profiles for the condensed and thermal clouds 
for typical values of gas parameters (ano = h/\/2miU). No = 
f d 3 rno(r) and Nt = / /riiT(r) are the total number of 
condensed and thermal atom respectively. 



while the non zero ones are determined from the semi- 
classical expression: 



1 



exp[/3(e^(r) + ^(r)- M )]-l 



(8) 



The set of Eqs. (|7l8| ) is reduced to a one dimensional 
problem if we assume the ansatz no(r) where r — 
^2mV{r) /uj and ZU = (w^WyWz) 1 / 3 . This ansatz is ex- 
act for a spherical trap and is accurate in the Thomas- 
Fermi limit Ui <C gn(0). It leads to the profiles in Fig|2] 
for the condensed and normal fluids and shows excellent 
agreements with both experiments and exact Monte- 
Carlo calculations in the determination of the den- 
sity profile of a trapped Bose condensed gas. These gen- 
eralizations allow the determination of the transferred 
momentum distribution N2 t k(t) = J d 3 r n2,k(r, t) from 
which we deduce the transferred thermal atom number 



EXPERIMENTAL PROPOSALS 
The gap experiment: four-photon process 

For q = and g±2 ~ g, the Raman spectrum becomes 
discrete in a homogeneous gas. The resonance frequen- 
cies correspond to a gap Swk,q(r) = —gn(r) associated 
to the exchange interaction energy for the single mode 
transition and to (gi2 — g)n(r) for the collective mode 
transition 



In comparison, if the condensed atom 



spectrum is quite similar, the thermal atom one displays 
differences in the Bogoliubov approximation. Since the 
energy difference £2,k(r) — e f k( r ) * s k dependant, no gap 
is observed and the oscillations are smoothed out leading 
to a continuous spectrum. 



In order to distinguish clearly between the discrete 
and the continuous spectra, the coupling potential acts 
specifically in the trap center in order to reduce inho- 
mogeneous broadening. This is realized by means of 
four beams (see FigJ3]) 18]: two gaussian astigmatic 
beams cr + polarized along the z axis of quantization 
with the intensity profile /i(r) = Ioi cxp(—2r 2 /wf(r z ) — 
2ry/w%(r z )) and two others ir polarized along the y axis 
with J 2 (r) = J 2 cxp(~2rl/wl(r y ) - 2r 2 z /wl{r y )) where 
Wi(s) = Wi(l + (sA) 2 /(tt 2 w;4)) 1 /2. The sum f tneir 

frequency differences corresponds to the transition fre- 
quency ui. Provided that A « uij, we define an effective 
waist w such that: 



U! W 



1.1 1 , 

— (^2+^2) 



LOyW^ 



LU Z W2 



(9) 



In these conditions, the resulting potential Vr(t) = 
Vfto exp(— 2r 2 /w 2 ) is optimized for an atom transfer in 
the most homogeneous region with q = 0. To fix the 
idea, we choose A = 843nm and w — 7/xm which reduces 
to about 10 4 the thermal atom effective number that can 
be specifically addressed. Transferring a small fraction 
of about 10% and for a detection resolution of about 100 
atoms, we obtain a signal to noise ratio of about 10. A 
relative difference in the scattering lengths is also needed 
to observe the gap resonance and is obtained from the ap- 
plication of an external magnetic field [HJ . These consid- 
eration leads to the spectra of FiglH Note the two orders 
of magnitude between the two peak intensities and the 
oscillatory behavior of period 1/t — lOQHz associated to 
the finite time resolution. The finite size of the beam 
provides an additional negligible frequency uncertainty 
of about hj (\/mf3w) in the resolution. 



The screening experiment: Raman scattering 

The absence of Raman transition due to screening is 
observed in the scattered atom momentum distribution. 
For a long time, the transient effects in Eq.@ can be 
neglected leading to a constant transfer rate and, except 
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FIG. 4: Transferred thermal and condensed atoms vs. the 
detuning in the GRPA and Bogoliubov approximation. 



for the fact that the external potential is screened, we 
recover the Fermi golden rule: 



^2,k+qM 



t 



^ 27rt3(r)n k (r)*(w 



- w k , q (r)) 



ft 2 |/C 12 (r,q, W )P 



(10) 



Considering 512 ~ g, the transition energy is position 
dependant causing inhomogeneous broadening: HSuj = 
k z q z /m + e q — gn(r). In the absence of screening, a reso- 
nance maximum appears for k z = 0. The screening factor 
strongly reduces the Raman scattering and forbids it at 

this maximum i.e. N2,k x ,k v ,q ie ('t)/t thus avoiding 
the condensed atom transfer. For simplicity, let uj x = uj y . 
The atoms are transferred by means of a Raman transi- 
tion resulting from two gaussian symmetric laser beams 
such that their wavevector difference q is along the z axis 
and their frequency difference is the transition frequency 
bj. For small q z , the angle between the beams is small 
and the Raman potential has the gaussian circular profile 
Vr = VRoexp(— 2(r x + ry)/w1{r z )). Once the atoms are 
transferred, the trap is switched off and after a time of 
flight, the density profile provides their momentum dis- 
tribution. 

A negative detuning is chosen in order to scatter the 
thermal atoms with k z positive in the trap center region 
and negative otherwise. The graphs in Figj5] illustrate 
well the total screening effect around k z = for which 
the macroscopic wave function deforms its shape in order 
to attenuate locally the Raman potential, thus preventing 
single atom scattering. The left part of the distribution 
(k z < — Q^im^ 1 ) shows the thermal atoms coming from 
the outer condensate region. The choice of q z is such that 
the LDA validity condition q z w$ ^> 1 is fulfilled but also 
such that, during the flight, the mean field energy does 
not affect much the momentum distribution. The inter- 
action time must be much lower than the relaxation time 
associated with collisions t <C t ~ \J (3m/8na 2 nT{0) to 
avoid the equilibrium relaxation of the momentum dis- 
tribution. Its finite value creates an energy uncertainty 
that alters the validity of Eq. (JTDJ by not suppressing to- 
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FIG. 5: Thermal atom distribution iV^fc, = X^a 



2,k ver- 



sus k z in presence of screening using Eq. (ll0[) (dashed line) and 
in absence (dotted line) and presence (full line) of screening 
taking into account finite interaction time corrections. 



tally atom scattering at k z = 0. Also, this time must 
be adequately chosen to suppress the condensed frac- 
tion due to the Rabi flopping associated to the collec- 
tive mode: Nz.att) = 2sin 2 [(Mw — £q)t/2] N 2 o,q where 



»2,q(t) 

A^ 20 , q ~ 2 / d 3 rn (r)[V R {r)/(h5u> 
the case of Fig[5J 



e q )] 2 = 0.19A^ 2 ,t for 



Phenomenological approach 

Although theoretical statements argue in favor of 
GRPA, we cannot exclude that none of the two approx- 
imations reproduces correctly the physical observation. 
In such a case, we can use a phenomenological approach 
assuming a transition process from an excitation of un- 
known energy e^ k to an excitation of energy e^k+q ano - 
a process of creation of two excitations of energy e^"_ k 



and 



2,k+q' 



Using a four photon process interacting in 



the uniform region of the gas, the fraction of scattered 
atoms is then written under the form analog to Eq.(|10|l: 



2 ,k+q z 
t 



± 



A ± (q Z7 k)S(u J ±ef^ + el M+cl ) (11) 



where A± (q z , k) represent the associated amplitude for 
such transition processes. Experimentally, the imag- 
ing in two dimensions allows only the determination of 
F(k x ,k z ) = J dk y N2.k+q^- Thus the quantity (fTTj) is 



determined from the Abel's transformation: 



N. 



2,k+q z 



dF(y,k z ) 
dy 



dy 



y 2 - (*2 + kl) 



=(12) 



By varying the parameters q z and w, the resonance po- 
sitions in the k space allow to reconstruct the dispersion 
relations e^ k and e^ k for the excitations. 
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CONCLUSIONS 

We explored the many body properties of a trapped 
Bose gas that can be extracted from a two-level hyper- 
fine transition in the GRPA and Bogoliubov approxima- 
tion. The calculated spectra not only show the existence 
of a second branch of excitation but also the total screen- 
ing of the external potential which prevents single con- 
densed atom transitions. If the external potential orig- 
inates from the presence of a thermal atom, this total 
screening prevents the binary collision between that ther- 
mal atom and any condensed one. In this scenario, the 
metastability of the relative motion between the normal 
and super fluids is explained by the absence of this ex- 
change collision process Q. The experimental observa- 
tion of these phenomena will improve our understanding 
of the exact nature of the elementary excitations and of 
the origin of metastable motions in superfluids. 
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